- 'AD-A033  989 


UNCLASSIFIED 


I OF  I ! 
AD  A „ 
033  989 


DAYTON  UNI V OHIO  RESEARCH  INST 
DAYTON  AIRCRAFT  CABIN  FIRE  MODEL. 
JUN  76  P M KAHUT 


O 


F/G  1/2 

VOLUME  III.  COMPUTER  PROGRAM  — ETC(U) 
DOT-F  A74WA-3532 
FAA-RD-76-120.3  NL 


A0AO33989 


010022 

Report  No.  FAA-RD-76-120,  I! 


DAYTON  AIRCRAFT  CABIN  FIRE  MODEL 


Volume  1 1 1 - Computer  Program  User's  Guide 


Peter  M.  Kahut 

University  of  Dayton 
Research  Institute 
Dayton,  Ohio  45469 


l 

! 

f 

i 

i 

t 

. 


June  1976 
Final  Report 


Document  is  available  to  the  U.S.  public  through 
the  National  Technical  Information  Service, 
Springfield,  Virginia  22161. 


D D C 
rjj[?f7Pr?:in  I7f? 

V 0£C  28  1976 

EEEnrrc; 

D 


. . 


Prepared  for 

U.S.  DEPARTMENT  OF  TRANSPORTATION 
FEDERAL  AVIATION  ADMINISTRATION 
Systems  Research  & Development  Service 
Washington,  D.C.  20590 

REPRODUCED  BY 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U.S.  DEPARTMENT  OF  COMMERCE 


* 


w 


This  document  is  disseminated  under  the  sponsorship  of 
the  Department  of  Transportation  in  the  iiterest  of  infor- 
mation exchange.  The  United  States  Government  assumes  no 
liability  for  its  contents  or  use  thereof. 


Technical  Report  Documentation  Page 


5.  Report  Dote 

June  1976 

6.  Performing  Organization  Coda 

8.  Performing  Organization  Report  No. 


10.  Work  Unit  No.  (TRAIS) 

11.  Contract  or  Grant  No. 

FA74WA-  3532 

13.  Type  of  Report  and  Period  Covered 

Final  Report 

July  1974  - March  1976 

14.  Sponsoring  Agency  Cooe 


1.  Report  No. 

FAA-RD-76-  120,  III 


j 4.  Title  and  Subtitle 

Dayton  Aircraft  Cabin  Fire  Model 

Volume  III  - Computer  Program  User's  Guide 

7.  Author's) 

Peter  M.  Kahut 


| 9.  Performing  Organization  Name  and  Address 

University  of  Dayton  Research  Institute 
300  College  Park  Avenue 
Dayton,  Ohio  45469 

j 12.  Sponsoring  Agency  Nome  and  Address 

j Department  of  Transportation 
j Federal  A viation  Administration 
j Systems  Research  and  Development  Service 
Washington.  D,  C.  20590  


| 15.  Supplementary  Notes 


16.  Abstract 

A basic  mathematical  model  and  computer  simulation  program  have  been 
developed  to  assess  the  smoke  and  toxic  gas  emissions  resulting  from  the 
burning  of  cabin  interior  materials  of  a wide- body  transport  aircraft  in  a 
full-scale  fire.  The  simulation  is  based  on  laboratory  test  data  on  the  cabin 
materials.  This  report  is  a guide  for  use  of  the  computer  simulation  program 
which  includes  instructions  for  input  data  preparation,  sample  input  and  output, 
basic  definitions  concerning  the  simulation  program  and  mathematical  model, 
and  a brief  description  of  the  program  structure.  This  report  consists  of 
three  volumes:  Volume  I is  entitled  "Basic  Mathematical  Model"  and 
Volume  II  is  entitled  "Laboratory  Test  Program". 

D D C 

r'<  DEC  29  1976 

jusisEinns 

D 


17.  Kay  Words  I 18.  Distribution  Statement 

aircraft  fire  safety,  aricraft  interior  Document  is  available  to  the  public 
materials,  smoke  and  toxic  gases,  air-  through  the  National  Technical  Infor- 
craft  cabin  fires,  wide- body  aircraft,  mation  Service,  Springfield,  Virginia 
room  fire,  mathematical  model,  22151 

comDuter  simulation,  fire  research 


Form  DOT  F 1700.7  (8-72) 


Reproduction  of  completed  poge  authorized 


19.  Security  Classif.  (of  this  report) 

20.  Security  Classif.  (of  this  page) 

21.  No.  of  Pages 

22.  Prie. 

UNCLASSIFIED 

UNCLASSIFIED 

52 

I 


METRIC  CONVERSION  FACTORS 


4 


This  contract  was  prepared  by  the  University  of  Dayton  Research 
Institute  for  the  Federal  Aviation  Administration  Systems  Research  and 
Development  Service  under  Contract  FA74WA-3532  during  the  period  July 
1974  to  March  1976.  The  report  describes  the  development  of  a mathe- 
matical model  of  a fire  within  the  cabin  of  a wide  body  commercial  transport 
category  aircraft.  The  report  is  divided  into  three  volumes  of  which  this 
is  the  third.  Volume  I,  entitled  "Basic  Mathematical  Model,  " describes  the 
development  and  presents  example  results  of  the  model.  Volume  II,  "Labor- 
atory Test  Program,  " presents  the  results  of  a laboratory  test  and  data 
collection  program  conducted  in  support  of  the  development  of  the  model. 
Volume  III,  "Computer  Program  User's  Guide,  " is  a guide  for  use  of  the 
computer  program  which  implements  the  mathematicaJ  model. 

This  contract  was  administered  under  the  direction  of  Mr.  Robert  C. 
McGuire  and  Mr.  Charles  C.  Troha  of  the  Systems  Research  and  Development 
Service,  ARD  520.  Work  was  performed  at  the  University  of  Dayton  under 
the  supervision  of  Mr.  Nicholas  A.  Engler,  supervisor  of  the  Applied  Systems 
Analysis  Division.  Other  personnel  at  the  University  who  have  contributed 
to  this  program  include  Mr.  James  K.  Luers,  Mr.  Jerry  B.  Reeves,  and 
Mr.  Charles  D.  MacArthur.  The  author  wishes  to  express  his  gratitude 
to  all  those  mentioned  for  their  support,  encouragement,  and  valuable 
technical  contributions.  The  author  also  wishes  to  thank  Ms.  Jacquelin 
Aldrich  and  Ms.  Peggy  Cummings  for  their  patient  assistance  in  preparing 
the  manuscript. 


TABLE  OF  CONTENTS 


SECTION  NO.  TITLE  PAGE  NO. 

1 INTRODUCTION  1 

2 BASIC  DEFINITIONS  3 

3 PROGRAM  CONSTRUCTION  9 

4 INPUT  DATA  PREPARATION  17 

4.  1 Input  Data  Cards  17 

4.2  Sample  Input  33 

5 PROGRAM  OUTPUT  39 

5.  1 Output  Options  39 

5.2  Output  Formats  39 

6 PROGRAM  STATISTICS  47 

7 PROGRAM  AVAILABILITY  49 


LIST  OF  FIGURES 

FIGURE  NO.  & TITLE 


2.  1 Location  of  Cabin  Lining  Surfaces 

2.2  Location  of  Seat  Groups 

2.  3 Seat  Groups  - Distinct  Surfaces 
4.  1 Numbering  of  Elements 


Preceding  page  blank 

vii 


J 


PAGE  NO. 


5 

6 
7 

31 


SECTION  1 


INTRODUCTION 

The  Dayton  Aircraft  Cabin  Fire  (DACFIR)  Simulation  Program  is  the 
computer  implementation  of  the  DACFIR  Model  developed  by  the  University 
of  Dayton  Research  Institute.  The  model  is  a set  of  equations  ana  logic 
designed  to  predict  the  time  history  of  the  build-up  of  smoke,  heat,  and  toxic 
gases  in  the  cabin  of  an  aircraft  subjected  to  fire  within  a period  representative 
of  po?t  crash  emergency  evacuation  time.  The  model  provides  a means  of 
tracking  the  development  of  the  fire  and  the  changes  in  the  cabin  envi  l onment 
with  time.  Input  to  the  model  includes  a description  of  the  cabin  geometry  and 
ventilation  conditions,  a description  of  the  material  properties  as  measured  by 
laboratory  tests,  and  a description  of  the  initial  fire  situation.  The  DACFIR 
Simulation  Program  as  presented  in  this  report  applies  to  a particular  cabin 
geometry,  that  of  a representative  wide- body  aircraft  cabin  section.  This  pro- 
gram has  been  designated  as  DACFIR,  Version  1,  May  1976.  A complete  de- 
scription of  the  DACFIR  ModeL  is  contained  in  Volume  I,  "Basic  Mathematical 
Model". 

The  intent  of  this  user's  guide  is  to  provide  instruction  for  the  efficient 
use  of  the  DACFIR  Simulation  Program.  This  documentation  contairs  basic 
definitions,  a flow  chart  and  description  of  the  main  program,  a description 
of  the  program  input  and  output,  sample  program  input  and  output,  program 
statistical  data,  and  information  concerning  the  availability  of  che  program 
code.  It  does  not  contain  a detailed  description  of  the  computer  program  code 
and  thus  is  not  intended  as  a complete  reference  source  for  the  computer 
programmer. 

*•  * I 

It  should  be  stressed  that  careful  preparation  of  the  input  cards  is  of 
the  utmost  importance.  A simulation  program  by  its  ^ery  nature  utilizes  a 
relatively  large  amount  of  computer  time,  and  seemingly  insignificant  errors 
in  the  input  can  easily  result  in  serious  errors  in  the  results  if  not  cause 
abnormal  termination  of  the  run. 
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SECTION  2 


BASIC  DEFINITIONS 


An  element  is  the  smallest  unit  of  surface  area  which  is  utilized  in 
the  simulation,  and  is  a square  whose  dimension  is  six  inches  per  side. 

An  element  may  exist,  at  any  specified  time,  in  one  of  four  primary  states, 
virgin,  smoldering,  flaming  or  charred  or  in  one  of  three  secondary  states 
which  represent  temporary  conditions  intermediate  to  the  primary  states. 

A surface  is  defined  as  consisting  of  a group  of  elements  all  of  which 
lie  in  the  same  horizontal  or  vertical  plane  and  whose  material  properties 
are  identical.  The  program  recognizes  twenty  (20)  cabin  lining  surfaces  and 
nine  (9)  seat  groups.  The  twenty  cabin  lining  surfaces  are  depicted  in 
Figure  2.  1 and  are  as  follows: 

1.  Carpet 

2.  Lower  Right  Sidewail  Panel 

3.  Right  Window  Reveals  and  Window  Transparencies 
(considered  one  surface) 

4.  Upper  Right  Sidewall  Panel 

5.  Right  Side  Passenger  Service  Unit 

6.  Right  Side  Stowage  Bin  Bottom 

7.  Right  Side  Stowage  Bin  Face 

8.  Right  Ceiling  Panel 

9.  Right  Center  Stowage  Bin  Face 

10.  Right  Center  Stowage  Bin  Bottom 

11.  Left  and  Right  Center  Passenger  Service  Units 

12.  Left  Center  Stowage  Bin  Bottom 

13.  Left  Center  Stowage  Bin  Face 

14.  Left  Ceiling  Panel 

15.  Left  Side  Stowage  Bin  Face 

16.  Left  Side  Stowage  Bin  Bottom 

17.  LeftSide  Passenger  Service  Unit 

18.  Upper  Left  Sidewall  Panel 

19.  Left  Window  Reveals  and  Window  Transparencies 
(considered  one  surface) 

20.  Lower  Left  Sidewall  Panel 

The  seat  groups  are  nine  (9)  in  number  and  are  referenced  as  shown 
in  Figure  2.  2. 
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1. 

1st  Row,  Left 

6. 

2nd 

Row, 

Right 

2. 

1st  Row,  Center 

7. 

3rd 

Row , 

Left 

3. 

1st  Row,  Right 

8. 

3rd 

Row, 

Cente  r 

4. 

2nd  Row,  Left 

9. 

3rd 

Row, 

Right 

5.  2nd  Row,  Center 


The  left  and  right  seat  groups  each  contain  three  individual  seats.  The 
center  seat  groups  each  contain  four  individual  seats. 

Each  seat  group  consists  of  seven  (7)  surfaces  as  shown  in  Figure  2.  3. 

1.  Cushion  Bottom 

2.  Lower  Rear  Backrest 

3.  Upper  Rear  Backrest 

4.  Backrest  Top 

5.  Backrest  Front 

6.  Cushion  Top 

7.  Cushion  Front 

The  computer  program  utilizes  seven  (7)  types  of  materials: 

1.  Carpet  Material 

2.  Sidewall  Material 

3.  Window  Reveal-Transparency  Material 

4.  PSU  Facing  Material 

5.  Stow  Bin  Material 

6.  Ceiling  Panel  Material 

7.  Seat  Upholstery  Material  with  Padding 

It  is  assumed  that  all  materials  can  yield  one  or  more  of  the  following 
toxic  gases  as  the  mat  ,rial  becomes  involved  in  the  fire: 


1. 

CO 

6. 

h2s 

2. 

HC1 

7. 

NH 

3. 

HCN 

8. 

NOx 

4. 

5. 

HF 

S°2 

9. 

COC1 

Whenever  surfaces,  seat  groups,  materials  or  toxic  gases  are 
referenced,  the  numbering  will  always  be  identical  with  that  given  above. 
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SECTION  3 

PROGRAM  CONSTRUCTION 


The  computer  simulation  program  was  written  in  modular  form. 

That  is,  each  subroutine  has  as  its  purpose  a specific  function  or  set  of 
computations.  These  subroutines  are  linked  by  means  of  the  coding  in 
the  'Main  Program'.  The  coding  in  the  'Main  Program'  then,  provides 
the  necessary  controls  for  the  logical  flow  of  the  sequence  of  computations. 

The  following  three  pages  contain  a flow  diagram  of  the  Main 
Program.  The  numbers  in  the  left  margin  are  utilized  to  reference  the 
diagram  to  a brief  description  of  each  program  step  or  block.  This 
narrative  is  contained  on  the  pages  immediately  following  the  flow  charts. 
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'MAIN  PROGRAM'  FLOW  CHART 


Component  No. 
wee  pp.  13-15) 


l. 


2. 


3. 

4. 

5. 


6. 


7. 

8. 


9. 


START 


« 

NPASS  = 0 
IPR1  = 0 
IPR2  = 2 

CALL  INPUTG 
(IN) 

(InI  0)? 25L*  PRINT: 

yes 

l 

PRINT:  'RUN  ID’ 


CALL  INPUTM 

I 

CA  LL  INIT 

I 

CALL  INPUTO 


A) ITIM2  = ITIM2  + ITSPRD 

1=0 
K = 0 
KPR  = 0 
ISW  = 0 
ISAVE  = 0 

ITIME  = ITIME  + IDELT 

I 

TDQ  = 0. 

TRSF  = 0. 

TRGF(J)  = 0;  J = 1,9 


'NORMAL 

COMPLETION 


—►STOP 
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'MAIN  PROGRAM’  FLOW  CHART  (Continued) 


Component  No. 
(see  pp.  13-15) 


1 6. a, b 


■1*1+  1 


NFEG)  - 0?- 


■ CALL  FIRE(I) 


K = KPR? 


CALL  RATES(I) 


I > LSN  ? 


CALL  COND(I) 
CALL  FCON(I) 
CALL  PVOL(I) 


KPR  = K 


CALL  CONDS I) 
CALL  FCONS(I) 
CALL  PVOLS(  I) 


O 


KPR  = K 


CALL  TEST  (I) 


■♦ISAVE  - 07- 
yes 

l 

■ ■-l  < NS? 


-© 


N FIRES  - K 


CALL  ELEM 


■ 


'MAIN  PROGRAM'  FLOW  CHART  (Continued) 


Component  No. 
(see  pp.  13-15) 


22. 

23. 

24. 


25. 


26. 

27 


28. 

29. 

30 

31. 

32. 

33. 

34. 


ISW 


CALL 


= 1 
AFP 


CA  LL  ATMOS 


NPASS  = NPASS+1 


yes 


NPASS  =1? 
no 

1 

IPR 1 = MOD(  ITIME,  IPEMS) 
IPR2  = MOD(  ITIME,  IPSPR) 


CALL  OUTPUT 


yes 


ITIME  > ITFIN?- 
no 

1 


ISW  = 19- 


yes 

l 

CALL  RESET 


ISW  = 0 


(ITIME  + IDELT)  > 1TIM2?- 
no 

1 

ITIME  = ITIME  + IDELT 


© 


PRINT: 


'NORMAL 

COMPLETION'' 


© 


STOP 
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COMMENTS  FOR  PROGRAM  FLOW  CHART 


No. 

1.  'NPASS'  is  a counter,  initialized  at  this  point,  which  will  contain 
the  number  of  times  the  cabin  atmosphere  computations  have  been 
performed.  IPR1,  IPR2  are  print  controls. 

2.  Subroutine  'INPUTG'  initializes  and  defines  those  variables  pertaining 
to  the  geometry  of  the  cabin  section. 

3.  Test  for  run  termination  (no  additional  input  cards  to  be  read). 

4.  Eighty  characters  of  run  identification  are  printed. 

5.  Subroutine  TNPUTM'  reads  all  input  data  pertaining  to  the  material 
properties  of  each  surface. 

6.  Subroutine  'INIT'  performs  basic  computations  irom  the  input  data 
and  initializes  those  variables  required  to  start  the  integration. 

7.  Subroutine  ’INPUTO'  reads  all  data  relating  to  the  ignition  source. 

8.  This  point  is  the  start  of  the  primary  integration  loop  in  the  program. 

TTIM2’  is  the  time  associated  with  the  flame  propagation  computations, 

’ITIME'  is  the  time  associated  with  the  cabin  atmosphere  computations. 

9.  Initialize  sums  which  contain  emission  rates  for  all  fires. 

10.  I'  is  the  surface  index  ( 1 through  20  are  cabin  lining  surfaces,  21 
through  29  are  seat  groups). 

11.  Test:  Any  flaming  elements  on  surface  T'?  If  not,  by-pass  flame 
propagation  computations  for  this  surface. 

12.  Subroutine  'FIRE'  isolates  a fire  on  the  specified  surface  and  performs 
the  computations  of  the  flame  properties. 

13.  If  K = KPR  at  this  point,  a new  fire  was  not  found  in  subroutine  'FIRE'. 

14.  Subroutine  ’RATES'  determines  the  heat  flux  at  various  points 

associated  with  one  specific  fire,  and  interpolates  for  material 
properties  as  a function  of  heat  flux. 

15.  If  I > LSN,  the  fire  under  consideration  is  located  on  a seat;  other- 
wise the  fire  is  located  on  a cabin  lining  surface, 

16.  Subroutines  'COND',  'CONDS'  determine  flame  propagation  via 
conduction. 
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COMMENTS  FOR  PROGRAM  FLOW  CHART  (Continued) 
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j 


I 


I 

■ 

: 

i 


No. 

Subroutines  'FCON1,  'FCONS'  determine  flame  propagation  via 
flame  c ontact. 

Subroutines  'PVOL1,  'PVOLS'  test  for  possible  elemental  change 
of  stace  due  to  the  pyrolysis  (smoldering)  of  elements  in  the  vicinity 
of  a fire. 

17.  Subroutine  'TEST'  determines  if  any  flaming  elements  change  to 
the  charred  state  and  sums  the  emission  rates  for  each  fire. 

18.  If  'ISAVE'  i 0;  then  return  control  to  subroutine  'FIRE'  (continue 
to  search  surface  ' I ’ for  fires). 

19.  Test:  Have  all  cabin  lining  surfaces  and  seat  groups  been  examined 
this  time  step?  If  not,  increment  for  next  surface. 

20.  The  variable  'NFIRES'  contains  the  total  number  of  distinct  fires 
in  progress  during  this  time  interval. 

21.  Subroutine  'ELEM'  updates  the  time  counters  and  indicators  associated 
with  each  element. 

22.  'ISW'  is  a switch.  When  ISW  = 1,  indicates  that  the  flame  propagation 
computations  have  been  performed  this  time  step. 

23.  Subroutine  'AFP'  determines  the  total  number  of  flaming  and  smoldering 
elements  and  sums  emission  rates. 

24.  Subroutine  'ATMOS'  contains  all  of  the  equations  describing  the  ca  :>in 
atmo  sphere. 

25.  Add  'one'  to  the  pass  counter. 

2b.  Test:  If  this  is  the  first  pass  through  the  program,  automatically 
print  flame  propagation  and  cabin  atmosphere  data. 

27.  Determine  if  flame  propagation  and/or  cabin  atmosphere  data  is  to  be 
printed  this  time  step. 

28.  Subroutine  OUTPUT'  consists  of  the  required  print  and  format 
statements  and  controls  to  obtain  the  output  data  as  required. 

29.  Test:  If  simulation  time  has  expired,  print  appropriate  message 
and  terminate  the  run. 

30. -  If  the  flame  propagations  computations  have  been  performed  this 

32.  pass,  reset  computer  words  containing  the  element  states  information. 
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COMMENTS  FOR  PROGRAM  FLOW  CHART  (Continued) 


No. 

33.  If  flame  propagation  computations  are  to  be  performed  the  next 
time  step,  re- enter  appropriate  Loop. 

34.  If  flame  propagation  are  not  required  next  time  step,  increment  cabin 
atmopshere  'At'  and  re-enter  cabin  atmosphere  computations. 
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SECTION  4 


INPUT  DATA  PREPARATION 

This  section  describes  the  input  requirements  of  the  DACFIR  Com- 
puter Program.  The  preparation  of  each  input  card  is  described,  and, 
where  necessary,  a brief  explanation  of  the  input  data  requirements  and 
options  is  included.  Familiarity  with  the  DACFIR  Mathematical  Model,  as 
described  in  Volume  I,  is  assumed.  Following  the  input  preparation 
instructions  is  a listing  of  a sample  input  data  deck.  The  specific  set  of 
input  data  shown  in  the  listing  was  used  to  create  Case  1 of  the  sample  runs 
discussed  in  Section  7,  Volume  1. 

In  the  data  description  shown  below,  three  format  types  are  refer- 
enced. They  are 

Type  Desc  ription 

A Alphanumeric,  any  combination  of  letters, 

numbers,  and  special  characters  (including 
blanks  ) may  be  entered  in  the  appropriate 
column. 

I Integer,  the  entry  must  be  right  justified 

in  the  field  (range  of  columns  ).  Example: 
when  the  number  '25'  is  entered  in  a five- 
column  field,  it  must  be  preceded  by  three 
blanks. 

F Floating  point,  the  entry  may  appear  any- 

where in  the  specified  field,  but  the  insertion 
of  a decimal  point  is  mandatory. 


4.  1 INPUT  DATA  CARDS 

Input  data  cards  are  to  be  prepared  as  described  below.  Some  of 
the  input  data  for  the  version  of  the  DACFIR  Program  described  here 
(Version  1,  30  May  1976)  is  to  be  regarded  as  fixed  and  not  user -defineable. 
This  data  consists  of  various  geometric  relationships  common  to  all  wide- 
body  cabin  interiors.  The  fixed  data  is  indicated  in  the  discussion  and  users 

Preceding  page  blank 
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of  the  program  should  prepare  this  data  according  to  the  example  data  set 
given  in  Section  4.  2.  Card  Number  in  the  list  refers  to  the  number  in  the 
sample  data  set. 


Card 

Card 

Format 

Type 

Numte  r 

Var. 

Dim. 

Col. 

Type 

Desc  ription 

I 

1 

IDENT 

20 

1 -80 

A 

Run  identification 

2 - 5 

2 - 133 

- 

- 

- 

- 

Fixed  input  data  (see 
sample  data,  pp.  33-  35) 

6 

134 

SQD 

- 

1-10 

F 

Element  square 
dimension  (fixed) 

6 

134 

RFWS 

* 

11-20 

F 

Flame  spread  rate 
sidewall  to  seat  or 
seat  to  sidewall  (ft/  sec) 

6 

134 

DWS 

■ 

21-30 

F 

Separation  distance 
outside  seat3  to 
sidewall  (ft) 

6 

134 

CH 

- 

31-40 

p 

Cabin  floor  to 
ceiling  height  (ft) 

6 

134 

CL 

■ 

41-50 

F 

Cabin  section  length 
(see  Figure  5.  1 f 
Volume  1 ) (ft ) 

6 

134 

C .V 

- 

51-60 

p 

Cabin  width  'ft) 

6 

134 

SL 

- 

61-70 

p 

Detailed  section 
length  (see  Figure  5.  1, 
Volume  1 ) (ft ) 

6 

134 

sw 

- 

71-80 

p 

Section  width  (mu3t  be 
equal  to  1 / 2 the  value  of 
C IV !( ft  Jin  columns  51-oG 

7 - 9 

135  - 138 

" 

- 

“ 

“ 

Fixed  input  (see  sample 
data,  p.  35) 

10 

139 

imatl 

20 

1-2 

39-40 

i 

Material  type  for  each 
surface  taken  in  order 
around  the  cabin  interior 
(material  type  denoted  by 
the  integers  1-7,  see  p.  4 

11 

140 

lmats 

7 

1-2 

13-14 

i 

Material  type  for  each 
surface  on  the  s »ats 
(see  Figure  2.  3 & p.  4) 

18 


Card 

Card 

Format 

Type 

Number 

Var. 

Dim. 

Col. 

Type 

De  serration 

12 

141 

GTAB 

7 

1-  10 

F 

Stochiometric  fuel  to 
oxygen  ratio  for  each 
combustion  of  tne  seven 

61-70 

material  types  (unitless) 

12 

142 

qtab 

7 

1-10 

F 

Effective  heat  of  com- 
bustion for  each  of  the 
seven  material  types 

61-70 

(Btu/lbm) 

12 

143 

rtab 

7 

1-10 

F 

Fuel  vapor  density  at 
base  of  fire  for  each 
of  the  seven  material 

61-70 

types  (lbm/ftl) 

12 

144 

utab 

7 

1-10 

p 

Fuel  vapor  flow  velocity 
at  base  of  fire  for  each 
of  the  seven  material 

61-70 

types  (ft/  sec  ) 

12 

145 

TP 

7 

1-10 

F 

Time  of  transition  of  an 
element  from  the  ambient 
to  the  smoldering  state, 

61-70 

each  material  type  (sec) 

12 

146 

TPC 

7 

1-10 

F 

Time  of  transition  of  an 
element  from  the  smol- 
dering to  the  charred 
3tate,  each  material 

61-70 

type  (sec ) 

12 

147 

RSS 

7 

1-10 

F 

The  smoke  production 
rate  for  each  of  the  seven 
material  types  in  the 
smoldering  state,  'par- 

61-70 

ticles'  / (ft-  ■ sec  ) 

(See  Volume  1,  Section  5 
for  discussion  of  the  units 
'particles  of  smoke'. 

12 

148 

RCS(  1 ) 

7 

1-10 

F 

Production  rate  of  CO  for 
each  material  in  the  smol- 
dering state,  microlbs/ 

61-70 

ft^  • sec  ) (1  microlb  = 10^ 

19 

Card 

Card 

Format 

Type 

Number 

Var. 

Dim. 

Col. 

Tvpe 

Description 

12 

U9 

RGS(2) 

7 

1-10 

F 

Production  rate  of  HC1 
for  each  material  in  the 
smoldering  state. 

61-70 

microlbs /(ft^  • sec) 

12 

150 

RGS(3) 

7 

1-10 

5* 

Production  rate  of  HCN 
for  each  material  in  the 
smoldering  state. 

61-70 

microlbs /(ft*2  • sec) 

12 

151 

RGS(4 ) 

7 

1-10 

F 

Production  rate  of  HF 
for  each  material  in  the 
smoldering  3tate, 

61  -70 

microlbs/(ft^  • sec) 

12 

152 

RGS( 5 ) 

7 

1-10 

F 

Production  rate  of  SOo 

for  each  material  in  the 
smoldering  state. 

61-70 

microlbs /(ft  ^ • sec) 

12 

153 

RGS(6) 

*7 

1-10 

TT 

Production  rate  of  H?S 
for  each  material  in  the 
smoldering  state. 

61-70 

microlbs / (ft^  • sec) 

12 

154 

RGS(7) 

7 

1-10 

F 

Production  rate  of  NH3 
for  each  material  in  the 

. 

smoldering  state. 

61-70 

microlbs/(ft“  • sec) 

12 

155 

RC-S(8 ) 

7 

1-10 

F 

Production  rate  Of  NO  for 
each  material  in  the 
smoldering  state, 

61-70 

microibs/(ft^  • sec) 

12 

156 

RGS(9) 

7 * 

1-10 

F 

Production  rate  of  CCCI— , 
for  each  material  in  the 

. 

smoldering  :tate. 

61-70 

microlbs/(ft“  ’ sec) 

Card  Type  13  (numbers  157  through  275  in  the  sample  deck)  consist 

of  tables  of  various  material  properties  as  a function  of  Q,  the  heat  flux  in 

2 

Btu/(ft  • sec).  Each  table  consists  of  six  X,  Y pairs  of  points  - - X being 
the  heat  flux,  Y the  material  property.  If,  during  the  simulation,  a material 
property  value  is  needed  for  a heat  flux  less  than  Q^,  then  the  value  Y^  will 


20 


be  used.  Linear  extrapolation  is  used  to  obtain  a material  property  value 
for  a heat  flux  greater  than  Q . The  formats  of  these  tables  are  as  follows. 


Cols. 

1-5 

6-13 

14-18 

19-26 

27-31 

32-39 

40-44 

Q1 

Y1 

Q2 

Y2 

Q3 

Y3 

Q4 

45-52 

53-57 

58-65 

66-70 

71-78 

Y4 

Q5 

Y 5 

Q6 

Y6 

It  is  not  necessary  that  the  heat  flux  values  be  separated  by  equal  incre- 
ments. All  values  are  input  in  'F'  format. 


Card 

1X22 

13 


Card 

No. 


157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 


Table  as  a Function  of  Heat  Flux  (Q) 

Horizontal  flame  spread  rate,  material  #1,  ft/sec 
Horizontal  flame  spread  rate,  material  #2,  ft/sec 
Horizontal  flame  spread  rate,  material  #3,  ft/sec 
Horizontal  flame  spread  rate,  material  #4,  ft/sec 
Horizontal  flame  spread  rate,  material  #5,  ft/sec 
Horizontal  flame  spread  rate,  material  #6,  ft/sec 
Horizontal  flame  spread  rate,  material  #7,  ft/sec 
Vertical  upward  flame  spread  rate,  material  #1,  ft/sec 
Vertical  upward  flame  spread  rate,  material  #2,  ft/sec 
Vertical  upward  flame  spread  rate,  material  #3,  ft/sec 
Vertical  upward  flame  spread  rate,  material  #4,  ft/sec 
Vertical  upward  flame  spread  rate,  material  #5,  ft/sec 
Vertical  upward  flame  spread  rate,  material  #6,  ft/sec 
Vertical  upward  flame  spread  rate,  material  #7,  ft/sec 
Vertical  downward  flame  spread  rate,  material  #1,  ft/sec 
Vertical  downward  flame  spread  rate,  material  #2,  ft/sec 
Vertical  downward  flame  spread  rate,  material  #3,  ft/ sec 
Vertical  downward  flame  spread  rate,  material  #4,  ft/sec 
Vertical  downward  flame  spread  rate,  material  #5,  ft/sec 
Vertical  downward  flame  spread  rate,  material  #6,  ft/sec 


I 
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Card 

Card 

Type 

No. 

Table  as  a Function  of  Heat  Flux  (Q) 

13 

177 

Vertical  downward  flame  spread  rate,  material  #7,  ft/sec 

1 1 

178 

Time  interval,  time  lag  from  the  time  of  flame  contact  to  the 
time  material  begins  flaming  combustion,  material  #1,  sec 

f 1 

179 

Time  interval,  time  lag  from  the  time  of  flame  contact  to  the 
time  material  begins  flaming  combustion,  material  #2,  sec 

11 

180 

Time  interval,  time  lag  from  the  time  of  flame  contact  to  the 
time  material  begins  flaming  combustion,  material  #3,  sec 

1 1 

181 

Time  interval,  time  lag  from  the  time  of  flame  contact  to  the 
time  material  begins  flaming  combustion,  material  #4,  sec 

If 

182 

Time  interval,  time  lag  from  the  time  of  flame  contact  to  the 
time  material  begins  flaming  combustion,  material  #5,  sec 

II 

183 

Time  interval,  time  lag  from  the  time  of  flame  contact  to  the 
time  material  begins  flaming  combustion,  material  #6,  sec 

11 

184 

Time  interval,  time  lag  from  the  time  of  flame  contact  to  the 
time  material  begins  flaming  combustion,  material  #7,  sec 

It 

185 

Heat  release  rate  while  material  #1  is  undergoing 
flaming  combustion,  BTU/(ft^  • sec) 

It 

186 

Heat  release  rate  while  material  #2  is  undergoing 
flaming  combustion,  BTU/(ft^  • sec) 

1 1 

187 

Heat  release  rate  while  material  #3  is  undergoing 
flaming  combustion,  BTU/(ft2  . sec) 

11 

188 

Heat  release  rate  while  material  #4  is  undergoing 
flaming  combustion,  BTU/(ft^  • sec) 

1 1 

189 

Heat  release  rate  while  material  #5  is  undergoing 
flaming  combustion,  BTU/(ft^  • sec) 

1 1 

190 

Heat  release  rate  while  material  # 6 is  undergoing 
flaming  combustion,  BTU/(ft2  . sec) 

rt 

191 

Heat  release  rate  while  material  #7  is  undergoing 
flaming  combustion,  BTU/(ft^  • sec) 

it 

192 

Smoke  production  rate  for  material  #1  in  the  flaming 
state,  'particles' /( ft^  • sec) 

(A  smoke  'particle'  is  the  amount  of  smoke  which  if 
contained  in  a volume  of  one  cubic  foot  would  cause  the 
light  transmission  over  a one  foot  path  to  be  reduced 
by  1 0%) 


22 


Card 

No. 


Table  as  a Function  of  Heat  Flux  (Q) 

Smoke  production  rate  for  material  #2  in  the  flaming 
state,  'particles' /(ff2  . sec) 

Smoke  production  rate  for  material  #3  in  the  flaming 
state,  'particles  /(ft^  . sec) 

Smoke  production  rate  for  material  #4  in  the  flaming 
state,  'particles  /( ft^  • sec) 

Smoke  production  rate  for  material  #5  in  the  flaming 
state,  particles' /(  ft^  . sec) 

Smoke  production  rate  for  material  #6  in  the  flaming 
state,  particles'  /( ft^  . sec) 

Smoke  production  rate  for  material  #7  in  the  flaming 
state,  'particles '/( ft2  . sec) 

Production  rate  of  CO  for  material  #1  in  the  flaming 
state,  microlbs/{  ft^  • sec) 

(1  microlb  = iG~®lb) 

Production  rate  of  CO  for  material  #2  in  the  flaming 
state,  mic  rolbs  /(  ft4*  • sec) 

Production  race  of  CO  for  material  #3  in  the  flaming 
state,  microlbs /( ft2  . sec) 

Production  rate  of  CO  for  material  #4  in  the  flaming 
state,  mic  rolbs  /( fm  • sec) 

Production  rate  of  CO  for  material  #5  in  the  flaming 
state,  microlbs /( ft2  . sec) 

Production  rate  of  CO  for  material  #6  in  the  flaming 
state,  mic  rolbs  /(  ft*-  • sec) 

Production  rate  of  CO  for  material  #7  in  the  flaming 
state,  mic  rolbs /( ft2  . sec) 

Production  rate  of  HC1  for  material  #1  in  the  flaming 
state,  mic  rolbs /(  ft^  • sec) 

Production  rate  of  HCl  for  material  #2  in  the  flaming 
state,  microibs  /( ft"  • sec) 

Production  rate  of  HCl  for  material  *3  in  the  flaming 
state,  microIbs/(ft2  • sec) 

Production  rate  of  HCl  for  material  #4  in  the  flaming 
state,  mic  rolbs /(  ft2  . sec) 

Production  rate  of  HCl  for  material  ^5  in  the  flaming 
state,  mic  rolbs /( ft2  . sec) 


I 


Table  as  a Function  of  Heat  Flux  ( Q) 


r 


Card 

Type 

Card 

No. 

Table  as  a Function  of  Heat  Flux  ( Q) 

13 

211 

Production  rate  of  HCl  for  material  #6  in  the  flaming 
state,  mic  rolbs /(  ft^  . sec) 

f 1 

212 

Production  rate  of  HCl  for  material  #7  in  the  flaming 
? ° 
state,  mic  rolbs /( ft“  • sec) 

t 1 

213 

Production  rate  of  HCN  for  material  #1  in  the  flaming 
state,  mic  rolbs /( ft^  . sec) 

1 I 

214 

Production  rate  of  HCN  for  material  #2  in  the  flaming 
state,  microlbs/fft^  . sec) 

11 

215 

Production  rate  of  HCN  for  material  #3  in  the  flaming 
state,  mic  rolbs /(  ft2  . sec) 

1 1 

216 

Production  rate  of  HCN  for  material  #4  in  the  flaming 
state,  microlbs/(  ft-  . sec) 

? 1 

217 

Production  rate  of  HCN  for  material  #5  in  the  flaming 
state,  mic  rolbs /( ft^  . sec) 

I 1 

218 

Production  rate  of  HCN  for  material  # 6 in  the  flaming 
state,  mic  rolbs /( ft^  • sec' 

1 1 

219 

Production  rate  of  HCN  for  material  #7  in  the  flaming 
state,  mic  rolbs /( ft-  • sec) 

1 1 

220 

Production  ra.e  of  HF  for  material  #1  in  the  flaming 
state,  mic  rolbs /(  fH  . sec) 

t f 

221 

Production  rate  of  HF  for  material  #2  in  the  flaming 
state,  mic  rolbs /( ft^  • sec) 

? t 

222 

Production  rate  of  HF  for  material  # 3 in  the  flaming 
state,  mic  rolbs  /(  fm  • sec) 

f 1 

223 

Production  rate  of  HF  for  material  ff 4 in  the  flaming 
state,  mic  rolbs  /(  ft-  • sec) 

1 1 

224 

Production  rate  of  HF  for  material  #5  in  the  flaming 
state,  microlbs/(ft  * sec) 

• 1 

225 

Production  rate  of  HF  for  material  #6  in  the  flaming 
state,  mic  rolbs  /( ft^  . sec) 

1 1 

226 

Production  rate  of  HF  for  ma  terial  #7  in  the  flaming 
state,  mic  rolbs  /( fm  • sec) 

1 1 

227 

Production  rate  of  SC>2  for  material  =fl  in  the  flaming 
state,  mic  rolos /( ft-  • sec) 

1 1 

228 

Production  rate  of  SOi  for  material  ?2  in  the  flaming 
state,  mic  rolbs /(  ft-  . sec) 

L 


Table  as  a Pune  ti  or  o i H ea t Flux  ( Q\ 

Production  rate  of  SO?  for  material  ff3  in  the  flaming 
state,  micro  Lbs  /(  ft*  . sec) 

Production  rate  of  SO?  for  material  #4  in  the  flaming 
state,  micro  lbs  /(ft*"  . sec) 

Production  rate  of  SO 2 for  material  #5  in  the  flaming 
state,  microl.bs/(l.t'-  • sec) 

Production  rate  of  SO?  for  materia.!  #6  in  the  flaming 
state,  microl.be /( .V-  . sec) 

Production  rate  of  SO?  for  material  #7  in  the  ilaming 
state,  microlbs /l ft^  • sec' 

Production  rate  of  H?S  for  material  rl  in  the  flaming 
state,  micro  Lbs  /( ft*’  • sec'1 

Production  rate  of  H2;S  for  material  $2.  in  the  flaming 
state,  microl.bs  /( ft“  • sec'i 

Production  rate  of  .T?S  for  material  #3  in  the  flaming 
state,  microlbs /(ft'*  . sec) 

Production  rate  of  .-I?S  for  material  #4  in  the  flaming 
state,  mi. c.  roibs  /( ft'-  . sec) 

Production  rate  of  H~,S  for  material  #5  in  the  flaming 
state,  micro  lbs /Of'  . sect 

Production  rate  of  P-.S  for  material  -to  in  the  flaming 
state,  micro  .os /(ft-  • sec) 

Production  rate  of  H?S  for  material  4~  in  the  flaming 
ita:e,  m icrolbn /(  :t  • . sec) 

P roduction  rate  o:  'MH3  for  material  #1  in  the  flaming 
state,  rmcrolbs/(  ft*  • sec) 

Production  rate  o'  \rH ? to-  material  #2  in  the  flaming 
, , , „ 3 

state,  mLcroibs/x  .t"  • sec) 

Production  rate  o : NHj  for  material  #3  In  the  flaming 
state,  micro  Lb  3 /(ft-  • sec) 

Production  rate  ct  NH3  for  material  #4  in  the  flaming 
state,  mi.crolbs/(  :t-  • sec) 

Production  rate  of  NH3  for  material  #5  in  the  flaming 
state,  micro  Lbs  /(  ft-  • sec) 


Production  rate  of  NH3  for  material  #6  in  the  flaming 
state,  mi c roib  3 /(  ft ' • sec) 


1 J 


fi 


Card  Card 

Ty>e  No. 

13  247 

" 248 

' ' 1. 49 

" 250 

" 251 

" 252 

" 253 

" 254 

" 255 

" 256 

" 257 

" 258 

" 2 59 

" 260 

" 261 

" 262 

" 263 


Table  as  a Function  of  Heat  F lux  (Q) 

Production  rate  of  NH3  for  material  #7  in  the  flaming 
state,  mic  rolbs /(ft^  • sec) 

Production  rate  of  NOx  for  material  #1  in  the  flaming 
state,  micro  lbs  /(ft4-  • sec) 

Production  rate  of  NGX  for  material  #2  in  the  flaming 
state,  micrclbs /( ft"  • sec) 

Production  rate  of  NOx  for  material  #3  in  the  flaming 
state,  microlbs /( ib*  • sec) 

Production  rate  of  NO  for  material  #4  in  the  flaming 
state,  microlbs /( ft*-  • sec) 

Production  rate  of  NOx  for  material  #5  in  the  flaming 
state,  microlbs /(ft*  • sec) 

Production  rate  of  NO...  for  material  #6  in  the  .Taming 
state,  microlbs/(  it-  • sec) 

Production  rate  of  NOx  for  material  #7  in  the  flaming 
state,  microlbs  /(  ft*  • sec) 

Production  rate  01  CGCl;?  for  material  #1  in  the  flaming 
state,  mic rolbs /( ft*  • sec) 

Production  rate  of  COCl?  for  material  #2  in  the  flaming 
state,  mic  rolbs  /(  ft*  * sec) 

Production  rate  of  COCl?  for  material  #3  in  the  flaming 
state,  mic  rolbs /( it"  » sec) 

Production  rate  of  COCL?  for  material  #4  in  the  flaming 
state,  microlbs  /(ft*  • sec) 

Production  rate  of  COCO  for  material  #5  in  the  flaming 
state,  mic  rolbs /(  fP'  • sec) 

Production  rate  of  COCl^  for  material  #6  in  the  flaming 
state,  microlbs/{  ft^  • sec) 

Production  rate  of  COCI2  for  material  #7  In  the  flaming 
state,  mic  rolbs  /(  ft*  • sec) 

Time  interval,  time  required  for  material  #1  to  stop 
smoldering  after  heat  flux  falls  below  the  threshold 
value  for  this  material,  sec 

Time  interval,  time  required  for  material  #2  to  stop 
smoldering  after  heat  flux  falls  below  the  threshold 
value  for  this  material,  sec 


.1 


26 


\ 


v.3rd  Card 

Type  No. 

12  254 

” 265 

" 266 

’ 267 

’ 2 58 

269 

’ 270 

" 271 

" 272 

273 

" 274 


Ts bie  a s ? F un ction  of  Kea t Flux  ( Q) 

Time  interval,  time  required  for  material  ^3  to  stop 
smoldering  after  heat  flux  falls  below  the  threshold 
value  for  this  material,  sec 

Time  interval,  time  required  for  material  #4  to  stop 
smoldering  after  heat  flux  falls  below  the  threshold 
value  for  this  material,  sec 

Time  interval,  time  required  for  material  **5  to  stop 
smoldering  after  heat  flux  falls  below  the  threshold 
value  for  this  material  sec 

Time  interval,  time  required  for  material  #6  to  stop 
smoldering  after  heat  flux  falls  below  the  threshold 
value  for  this  material,  sec 

Time  interval,  time  required  for  material  fr7  to  stop 
smoldering  after  heat  Hue  falls  below  the  threshold 
value  for  this  material  : ec 

Time  interval,  from  time  material  #1  begins  flaming 
combustion  to  time  material  becomes  charred,  sec 

Time  interval,  from  time  material  #2  begins  flaming 
combustion  to  time  ma  erial  becomes  charred,  sec 

Time  interval,  from  time  material  #3  begins  flaming 
combustion  to  time  material  becomes  charred,  sec 

Time  interval,  from  time  material  4 4 begins  flaming 
combustion  to  time  mate;  iai  becomes  charred,  sec 

Time  interval,  from  time  material  45  begins  flaming 
combustion  to  time  material  becomes  charred,  sec 

Time  interval,  from  time  material  #6  begins  flaming 
combustion  to  time  male:  iai  becomes  charred,  sec 

Time  interval,  from  time  material  4 7 begins  flaming 
combustion  to  time,  material  becomes  charred,  sec 


n — 

- t 


Two  methods  may  be  used  to  specify  the  ignition  source  for  a simu- 
lation. The  first  method  is  illustrated  by  the  sample  data  of  Section  4.  2 and 
the  data  preparation  instructions  are  giver,  immediately  below.  This  first 
method  consists  of  specifying  the  amount  and  location  of  an  ignition  source 
material,  such  as  a flammable  liquid.  The  ignition  source  material  must 
be  located  on  one  of  the  cabin  interior  surfaces  and  must  be  described  by 
the  flammability  and  combustion  toxicity  parameters  given  below.  At  the 
start  of  the  simulation,  all  of  the  ignition  source  material  is  assumed  to  be 
on  fire  and  this  fire  will  continue  to  burn  for  a period  of  time  determined 

by  the  expression  t.  = m,. /m,  = m „./{p  .u  .Ar.)  where  t.  is  the  burn  time 

1 e lg  fi  ft  II  oi  oi  ft  lg 

in  seconds,  m...  is  the  mass  of  ignition  source  material  in  lbm,  m is  the 
fi  ' tt 

material  burning  rate  which  is  the  product  of  p , the  material  fuel  vapor 
density  entering  the  flame  base,  u the  fuel  vapor  velocity  entering  the 
flame  base,  aid  A^.  the  total  area  of  the  ignition  source  material.  For 
further  discussion  of  this  method  cf  describing  ignition,  see  Section  7 of 
Volume  I.  The  second  method  of  specifying  ignition  is  discussed  on  page  32 
of  this  volume. 


Card 

Type 

14 

Card 

No. 

Var. 

Sim. 

Col. 

T vne 

Definition 

276 

QP 

7 

i-  10 
61-70 

F 

Heat  flux  ( threshold  value)  at 
which  each  material  type  begins 
to  smolder,  Btu/(ft^  • sec) 

15 

277 

C PM 

-- 

1-  10 

F 

Specific  heat  oi  materials  at 
ambient  condi. ions,  average 
value,  Btu/(lbm  • °R) 

RHOM 

-- 

21-30 

F 

Bulk  density  of  materials, 
average  value  (lbm/ft^) 

XK 

-- 

21-30 

jp 

. Thermal  conductivity  of 
materials  at  ambient  conditions, 
average  value,  Stu/(ft  • sec  * ° 

XPEN 

-- 

31-40 

F 

Heat  penetration  depth  of 
materials,  average  value,  ft 

TO 

_ _ 

41-50 

p 

Ambient  temperature,  °R 

28 


Col. 


rehnition 


I 


l 


ij 


Card  Card 

Type  No.  Var.  Dim. 

16  278  DELTAT  --  1-10  F Integration  t ime  interval  for 

the  cabin  atmosphere  computa- 
t ons,  sec 

Tl'INAL  --  11-2  2 F End  time  lor  the  simulation 

run,  sec 

IRATXO  --  21-23  Ratio:  integration  time  interval 

of  flame  spread  computations 
integration  time  interval  of 
cabin  atmosphere  computations 
'must  be  > 1) 


IPEMS 

26-30 

* 

Print  interval,  caoin  atmosphere, 
sec 

IPSPR 

31-36 

r 

Print  interval,  flame  spread 
data,  sec 

17 

279 

NV 

1-5 

T 

Number  of  open  doorways 
{ may  be  = 0) 

BVENT  -- 

12-20 

s' 

Distance,  floor  to  top  of  open 
doorway,  it  (same  for  all 
doorways) 

18 

280 

WVENT  -- 

* 

1 • 1 

*/»  ^ 

o o 

Width  of  each  open  doorway,  ft 
*This  card  must  be  deleted 
if  NV  = 0 

19 

281 

QBKGND  -- 

l-iO 

y 

Background  heat  flux.  Btu/(ft2<  sec) 

20 

282 

IGSN 

1-5 

i 

Surface  on  which  ignition  source 
material  resides.  If  this  value  is 
entered  as  zero,  the  second  method 
of  specifying  ignition  is  indicated, 
ic  this  case,  no  cards  of  type  21,  22, 
or  23  should  appear  in  this  data  deck. 

QCI 

6-15 

- 

miective  heat  of  combustion 
of  the  ignition  source  material, 
£tu/lbm 

GAMI 

lb-25 

F 

Stochiometric  fuel  to  oxygen 
ratio  of  ignition  source  material 

RHOI 

26-  35 

F 

Fuel  vapor  density  at  base  of 
ignition  source,  lbm/ft-^ 

XMUI 

36-45 

r 

Fuel  vapor  flow  velocity  at 
base  of  ignition  source,  ft/sec 

XMFI 

46-55 

F 

Mass  of  ignition  source 
material  to  be  burned,  Ibm 

Card 

Card 

Type 

No. 

Var. 

Dim. 

Col. 

Type 

F 

Description 

21 

283 

RSI 

— 

1-10 

Smoke  production  rate,  ignition 
source  material,  ’pa  rticl»3!  /ft^  • sec 

RTGI(l) 

-- 

11-20 

F 

production  rate  of  CO,  ignition 

source  material,  microlb/(ft^  • sec) 
(1  microlb  = 10“^  Ibm) 

RTGI(2) 

-- 

21-30 

F 

production  rate  of  HCl,  ignition 
source  material,  microlb/(ft2  . sec) 

RTGI(3) 

-- 

31-40 

F 

Production  rate  of  HCN,  ignition 
source  material,  microlb/tft^  . sec) 

RTGI(4) 

41-50 

F 

Production  rate  of  HF,  ignition 
source  material,  microlb/(ft2  . sec) 

RTGI(5) 

-- 

51-60 

F 

Production  rate  of  SO2,  ignition 
source  material,  microlb/(ftr  • sec) 

RTGI(6) 

-- 

61-70 

F 

Production  rate  of  HgS,  ignition 
source  material,  microlb/fft^  . sec) 

RTGI(7) 

*- 

71-80 

F 

Production  rate  of  NHj,  ignition 
source  material,  microlb/(ft^  . sec) 

21 

284 

RTGI(8) 

-- 

1-10 

F 

Production  rate  of  NO^,  ignition 
source  material,  micfotb/(ft^  • sec) 

RTGI19) 

-- 

11-20 

F 

Production  rate  of  COClj,  ignition 
source  material,  microlb/(ft^  • sec) 

22 

285 

NIJSQ 

-- 

1-5 

I 

The  number  of  elements  covered 
by  the  ignition  source 

PIGN 

6-15 

p 

Perimeter  of  ignition  source 
fire,  ft 

23 

286 

IGNIJ  ( 

2100)  | 

j 1-5 

I 

i index  of  ignition  source  element 

thru 

301 

1 

1 6-10 

I 

j index  of  ignition  source  element 

Enter  one  pair  of  i,  j indices  per 
card  (total  number  of  cards  will  be 

equal  to  the  value  of  NIJSQ) 
For  numbering  of  elements, 
see  Figure  4.  1 . 

24 

302 

NIJC 

1-5 

I 

Number  of  elements  (on  any  surface  1 
to  be  set  to  the  charred  (inert)  state 
at  the  start  of  the  simulation.  If  the 
value  = zero,  do  not  include  cards  of 
type  25  in  the  deck. 

25 

(not 

I,  J 

1 

(1-5 

i index  of  inert  element 

shown ) 

< 

( 6-10 

j index  of  inert  element 

Enter  one  pair  of  i,  j indices  per 
card  (total  number  of  cards  will  be 
equal  to  the  value  of  NIJC 


] 
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The  second  method  for  specifying  the  ignition  source  is  somewhat 
more  simple  than  the  first.  In  this  method,  any  number  of  elements  on  one 
or  more  surfaces  can  be  set  to  the  flaming  state  at  the  start  of  the  simulation. 
To  indicate  this  method,  a zero  is  entered  in  column  5 of  card  type  20. 
Initialization  of  the  elements  to  the  flaming  state  is  done  by  entering  values 
on  card  types  26  and  27. 


Card 
T /pe 

Card 

No. 

Var.  Dim. 

Col. 

Tvpe 

Cesc  ription 

26 

303 

NOFL  -- 

1-5 

I 

•Total  number  of  flaming  elements 
on  surface  ISFL. 

ISFL 

6-10 

Surface  number  on  which  the  elements 
identified  on  the  following  type  27 
cards  are  located. 

27 

(not 

I.  J 

jl-5 

I 

i index  of  flaming  element 

shown) 

(6-10 

j index  of  flaming  element 

Enter  one  pair  of  i,  j indices  per 
card  to  a total  of  NOFL  cards. 

If  elements  are  to  be  initialized  to  flaming  on  additional  surfaces,  the  deck 
should  continue  with  groups  of  cards  of  types  26  and  27  arranged  as  above 
until  all  desired  elements  are  initialized.  No  further  cards  are  necessary 
following  the  last  type  27  card. 


*If  the  value  of  NOFL  entered  is  zero,  a zero  must  be  entered  for  the 
variable  ISFL  and  this  card  will  be  the  last  card  in  the  deck.  No  initiali- 
zation to  the  flaming  state  will  be  performed. 


J 
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4.2  SAMPLE  INPUT 


s 


The  following  six  pages  present  a sample  input  data  deck. 


o 

U 

Card  Card 

T ype  No. 

1 1 WI0E-8O0Y  CASE  1 S 0 A MATERIALS,  FLOOR  SPILL  EIRE 


2 

20 

<3  7 

22 

3 

0. 

0, 

1.0 

II 

1.0 

0. 

3 . 

II 

5 

1.0 

0. 

0 . 

It 

1.0 

0. 

0. 

II 

0. 

0. 

-1.0 

II 

0. 

0. 

-1.0 

II 

l.G 

0. 

0 . 

II 

10 

0. 

0. 

-1.0 

II 

-1.0 

0. 

0 . 

II 

0. 

0. 

-1.0 

II 

n. 

0. 

-1,0 

It 

o . 

0. 

-1.0 

It 

15 

1.0 

0. 

0 . 

II 

0 . 

0. 

-1.0 

It 

-1.0 

0. 

0. 

II 

0 . 

0. 

-1.0 

II 

0. 

0. 

•1.0 

II 

20 

-1.0 

0. 

0. 

It 

-1.0 

0. 

0. 

3 

-1.0 

0. 

0. 

4 

1 

0 1 

A0  1 

15 

0 .0 

5 

21 

1 10 

2 5 

1 

0 

1 

-1 

0 

5 

25 

22 

1A  27 

1 A 

1 

-13 

1 

0 

0 

5 

23 

31  A0 

2 5 

1 

-30 

1 

-1 

0 

5 

2>* 

1 10 

7 13 

1 

0 

1 

-6 

0 

5 

25 

1A  27 

6 O 

1 

-13 

1 

-5 

3 

5 

26 

31  A0 

7 10 

1 

-30 

1 

-6 

0 

5 

30 

27 

1 10 

12  15 

1 

0 

1 

-11 

0 

5 

23 

1 A ?7 

11  1A 

1 

-13 

1 

-10 

0 

5 

23 

31  A0 

12  15 

1 

-30 

l 

-11 

0 

4 

9 

A A 1 

A3  1 

15 

0.0 

5 

5 

A 1 5 A 

1 15 

0 

53 

1 

0 

0 

5 

35 

-23 

A3  A3 

2 A 

0 

10 

-1 

23 

0 

5 

-26 

A3  A3 

7 0 

0 

10 

-l 

28 

0 

5 

-29 

A3  A3 

12  1A 

0 

10 

-1 

33 

0 

4 

0 

7 AA 

51  1 

15 

0.0 

5 

5 

A 1 5A 

1 15 

0 

55 

1 

0 

0 

5 

40 

-23 

AA  AO 

5 5 

0 

10 

-1 

62 

1 

5 

-23 

AA  A0 

5 5 

0 

10 

1 

-38 

1 

5 

-26 

AA  A 9 

13  10 

0 

10 

-1 

62 

1 

5 

-26 

AA  A9 

10  10 

0 

IQ 

1 

-33 

1 

5 

-20 

AA  A 9 

15  15 

0 

10 

-1 

62 

1 

5 

45 

-20 

AA  A0 

15  15 

0 

IQ 

1 

-33 

1 

4 

0 

1 52 

5A  1 

15 

0.0 

33 


lT> 

0 

U 

Card 

Card 

Type 

No. 

41 

54 

1 

15 

55 

5 

47 

5 

0 

1 

0 

0 

4 

0 

0 

55 

56 

1 

15 

7.0 

4 

0 

0 

57 

60 

1 

15 

7.0 

4 

50 

1 

1 

61 

62 

1 

15 

0.0 

c 

a 

61 

62 

1 

15 

C 

63 

1 

0 

a 

D 

A 

0 

0 

63 

71 

1 

15 

8.0 

4 

0 

1 

72 

73 

1 

15 

0.0 

5 

4 

8 

72 

73 

1 

15 

C 

71 

1 

0 

a 

55 

0 

0 

74 

77 

1 

15 

7.0 

4 

0 

0 

73 

79 

1 

15 

7.0 

A 

0 

0 

80 

93 

1 

15 

7.0 

*X 

A 

0 

1 

34 

85 

4. 

15 

0.0 

t 

5 

4 

14 

64 

35 

1 

15 

0 

86 

1 

a 

0 

60 

0 

0 

86 

94 

1 

15 

8.0 

4 

0 

1 

95 

96 

1 

15 

Q.C 

5 

14 

96 

96 

1 

15 

a 

94 

1 

a 

0 

4 

0 

a 

97 

100 

1 

15 

7.0 

4 

0 

0 

101 

102 

1 

15 

7.0 

4 

65 

0 

l 

103 

105 

1 

15 

0.0 

5 

17 

103 

116 

1 

15 

C 

102 

1 

a 

0 

4 

o 

7 

106 

113 

1 

15 

0.0 

5 

17 

103 

116 

1 

15 

G 

102 

1 

0 

a 

5 

-21 

108 

113 

5 

5 

0 

1 

1 

-ss 

l 

5 

70 

-21 

108 

113 

5 

5 

0 

1 

-1 

119 

l 

5 

-24 

106 

113 

10 

10 

C 

1 

1 

-95 

i 

5 

5 

-24 

108 

113 

10 

1C 

0 

4 

A 

-1 

119 

i 

-27 

103 

113 

15 

15 

0 

1 

1 

-95 

l 

5 

-27 

106 

113 

15 

15 

0 

1 

-1 

119 

l 

4 

75 

0 

4 

114 

116 

1 

15 

0.0 

5 

17 

103 

116 

1 

15 

0 

10  2 

1 

0 

0 

5 

-21 

114 

114 

2 

4 

G 

1 

-1 

23 

3 

5 

-24 

114 

114 

7 

9 

0 

1 

-1 

28 

0 

5 

-27 

114 

1 14 

12 

14 

0 

1 

-1 

33 

a 

4 

80 

1 

6 

1 

10 

4 

22 

1.0 

5 

16 

3 

6 

5 

18 

-1 

10  3 

0 

5 

0 

5 

17 

1 

2 

5 

18 

-1 

10  3 

0 

5 

0 

5 

-19 

1 

1 

5 

11 

-1 

119 

0 

5 

i 

5 

-19 

1 

1 

13 

18 

1 

95 

0 

4 

i 

5 

85 

-20 

1 

1 

1 

4 

C 

114 

1 

1 

0 

5 

-20 

1 

1 

19 

21 

0 

114 

-l 

23 

a 

4 

1 

3 

1 

14 

1 

22 

1.0 

5 

12 

3 

6 

5 

16 

-1 

86 

a 

4 

0 

5 

11 

7 

3 

5 

18 

-1 

66 

0 

4 

0 

5 

90 

10 

9 

12 

5 

16 

-1 

86 

0 

4 

G 

4 

1 

6 

1 

10 

1 

22 

1.0 

5 

5 

9 

10 

5 

13 

-1 

65 

0 

5 

0 

5 

6 

5 

6 

5 

If 

-1 

65 

0 

5 

0 

5 

-3 

10 

10 

5 

11 

1 

38 

0 

5 

1 

5 

95 

-3 

10 

10 

13 

16 

-1 

62 

0 

4 

1 

5 

-2 

10 

10 

1 

4 

* 

V 

43 

1 

1 

0 

5 

-2 

10 

10 

19 

21 

Q 

43 

-1 

23 

0 

4 

1 

6 

1 

10 

1 

22 

1 • w 

1 


34 


L 


— 


o 

u 

Card  Card 


Type  No. 


10 

l 2 

3 2 4 5 5 6 

5 5 4 5 

5 6 5 5 4 2 3 

2 

11 

140 

7 7 

7 7 7 7 7 

12 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

12 

7 J 30 

. 7300. 

7000. 

7430  . 

730 

3 • 

7000  . 

7000 

12 

C < 0935  0.0  935 

0.C935  0.0935 

0 .0' 

935 

0.0935 

3 .0935 

12 

0.0374  0.0374 

0.0374  3.0374 

0.0374 

0 • 3374 

0.0374 

145 

9. 

20. 

12. 

13. 

14. 

7. 

5. 

16 

17 

500. 

67. 

150. 

15C. 

130 

• 

120. 

6QG. 

16 

12 

13. 

25. 

25. 

25. 

25.- 

75. 

10  . 

12 

2 • 

33. 

16. 

16. 

30  . 

4. 

10  . 

12 

0.7 

4. 

0.1 

0.2 

4. 

9 . 74 

C.  84 

12 

150 

0.32 

0.05 

0. 

0 . 

0.05 

0. 

0.62 

12 

0. 

6. 

0. 

3 . 

6. 

2. 

3. 

12 

3.0 

0. 

0. 

0 . 

0 . 

0 . 

3. 

12 

0. 

0. 

0 . 

0 . 

0. 

3. 

0. 

12 

0. 

0. 

0. 

0 . 

0 . 

0. 

0. 

12 

155 

0. 

0. 

0 . 

0. 

0 . 

0. 

3. 

12 

0. 

0. 

G. 

0. 

0. 

3 « 

0. 

13 

.3 

3.  1.23 

.017 

1 .54  .0  46  4 

2.32 

. 139 

3,96 

.417 

4.3 

.5 

n 

0. 

3.  1. 

0. 

2.2  . 0 C 9 3 

3.  06 

.0161 

4.0 

. 0 275 

5.0 

.038 

it 

1. 

0.  2.2 

.002 

3.06  • OC  4 2 

4.41 

. OOS 3 

4.5 

.0094 

3.0 

.312 

ft 

160 

1.32 

0.  2.2 

.0016 

3.05  .0034 

4.41 

.3075 

4 • 6 

. 0 086 

5.3 

.Gil 

m 

0 . 

0.  1.32 

0. 

2 . 1 0 .0155 

3.0  8 

.0131 

4.0 

. 0275 

3.0 

.038 

it 

0. 

3.  1. 

0. 

1.9  .0106 

2.82 

.0213 

3.96 

.3415 

4.20 

.045 

M 

0.5 

0.  1.23 

.004 

1.94  .0115 

3 . 5Q 

. 0660 

4 . 0 

.097 

5.0 

. 1 7 

it 

0 . 

0.  1.0 

0. 

2.0  0. 

3. 

3. 

4, 

C. 

5 . 

0. 

ti 

165 

5.4 

0.  1. 

0. 

2.2  .0192 

3.0  8 

.03 

3.5 

. G 35 

4 , 

.041 

ii 

0. 

0 . 1 . 

0. 

2.  .007 

3.06 

.0065 

4 , 4 i 

.0133 

5. 

.3175 

ii 

0. 

0.  1.32 

0 . 

2.2  .0035 

3.03 

. 0044 

4.41 

.0113 

5. 

.015 

ii 

0. 

3.  1. 

0. 

2.2  O'. 0152 

3.  06 

0.030 

3.5 

0.335 

4.0 

0.041 

ii 

0. 

0.  1. 

3. 

2.  0. 

3 • 

u , 

4. 

0 , 

5. 

0. 

ti 

170 

0.5 

0.  1.23 

C .006 

1.94  .0173 

3. 

.045 

3.96 

.08  34 

5. 

.15 

n 

0. 

3.  1. 

c . 

2.  0. 

3. 

0 . 

4. 

0 « 

5. 

0 . 

n 

0. 

0.  1. 

c. 

2.2  .00775 

7.03 

.0167 

3.5 

.0  24 

4 . 

. 0 305 

it 

0 . 

0.  1. 

0. 

2.2  .0116 

3.0  8 

. 0038 

4 . *•  1 

. 0 3 83 

5. 

.3311 

tt 

3. 

C.  1.32 

c. 

2.2  .2.14 

3.08 

. 0 0 34 

4.41 

. 0 3 75 

5. 

.0097 

it 

175 

0.0 

3.  1. 

0. 

2.2  0,3079 

3.38 

3.0197 

3.5 

0.C24 

4.0 

0.0305 

ii 

0. 

0.  1. 

0. 

2 . u . 

3. 

0. 

4. 

0. 

5. 

3. 

ii 

0.5 

9.  1.23 

. 33  32 

1.94  .Cl  9 2 

3 , 

.3  24 

3.36 

• 0 **4  5 

5. 

.03 

ii 

0. 

4000.  1. 

10. 

2.  5. 

3. 

2. 

4. 

1. 

5C. 

3. 

it 

0. 

9000.  0.5 

5 • 

2.2  1.92 

3.08 

4 , 44 

6. 

0.5 

50. 

0 . 

tt 

180 

0 . 

9300.  1. 

5C. 

2.2  21.6 

3 . 0 o 

10. d 

4.41 

2.43 

53. 

0. 

it 

0 . 

9303.  1. 

50. 

2.2  21.5 

3.33 

10  .8 

4.4; 

2.43 

5‘j . 

G. 

it 

c. 

9000.  1. 

10. 

2.  5. 

3 . 

2. 

4. 

A . 

50. 

0. 

it 

c. 

9000.  0.5 

5 . 

2.2  1.12 

3.08 

1.44 

6. 

0.5 

50 , 

0. 

tt 

0 . 

90  OC . 1. 

10 . 

2.2  4. 

4, 

1 , 

8. 

0.5 

50. 

0. 

tt 

185 

0. 

3 . 0.4 

0. 

1. 94  3.73 

2.82 

7. 06 

5. 

10  . 

6. 

10.4 

tt 

c. 

3.  3.5 

0. 

2.2  2.7  5 
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5.  1 OUTPUT  OPTIONS 

The  user  has  the  capability  of  specifying  the  frequency  of  the  output 
for  both  flic  cabin  atmosphere  and  the  flame  spread  data.  Input  card  type 
16  contains  the  required  parameters,  IPEMS  and  IPSPR.  'IPEMS',  the 
print  time  interval  for  the  cabin  atmosphere  data,  is  ordinarily  set  to 
correspond  to  every  tenth  pass  through  the  program  therefore,  if  the  inte- 
gration time  step  is  one  second,  set  'IPEMS'  = 10  seconds.  'IPSPR'  trie 
print  time  interval  for  the  flame  spread  data,  is  ordinarily  set  to  either 
'IPEMS'  or  3X  IPEMS.  The  maximum  flame  spread  rate  utilized  in  the 
computer  run  should  be  considered  in  determining  a value  for  'IPSPR'. 

Note  - Both  outputs  will  always  be  printed  upon  the  completion  of  the  first 
pass  through  the  program,  regardless  of  the  values  of  'IPEMS'  and  'IPSPR'. 


5.2  OUTPUT  FORMATS 


1.  The  following  will  appear  in  the  upper  left  corner  of  each  page 

of  output. 


TIME  = XXXX  SEC  AFTER  IGNITION 


■mere  XXXX  is  the  number  of  seconds. 

2.  Cabin  Atmosphere  Summary  (one  page  of  output) 

The  cabin  atmosphere  summary  consists  of  zone  deptn, 
gas  density,  gas  temperature,  material  surface  temper 
ature,  and  the  heat  rate  to  the  surface.  These  variables 
are  printed  for  both  the  upper  and  lower  zones. 

The  flow  rate  (in  and  out)  through  open  doorways  is  pr<  - 
sented  along  with  the  upward  gas  flow  resulting  from  all 
fires  in  the  cabin  section. 

The  smoke  concentration  and  toxic  gas  concentrations 
(for  nine  toxic  gases)  are  presented  in  two  forms:  the 
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values  computed  assuming  that  all  smoke  and  toxic  gases 
are  contained  within  the  upper  gas  layer,  "stratified", 
and  values  computed  assuming  that  the  smoke  and  toxic 
gases  are  uniformly  mixed  over  the  entire  cabin  section 
volume,  "uniform  mixing". 

The  values  reported  for  the  stratified  condition  should  be 
regarded  as  the  best  estimate  of  the  actual  concentrations 
achieved.  These  values  are  computed  by  the  methods 
described  in  Volume  I.  The  uniform  mixing  values  are 
computed  within  the  output  subroutine  by  multiplying  the 
stratified  values  by  the  ratio  of  the  upper  zone  depth  to  the 
cabin  height.  These  values  represent  a possible  lower 
limit  on  the  concentrations  which  would  result  if  the  air- 
borne materials  are  uniformly  distributed  over  the  entire 
available  volume. 

Two  columns  of  information  to  the  right  of  each  toxic  gas 
concentration  indicate  whether  the  concentration  has 
exceeded  a present  short  exposure  irritation  (IR)  level  or 
short  exposure  life  danger  (LD)  level.  A value  of  1 in  the 
appropriate  column  indicates  that  the  level  has  been  exceeded. 
The  preset  levels  used  are 


Irritation 

Life  Danger 

(ppm) 

(ppm) 

CO 

1.  0 

10,  000. 

HC1 

35.  0 

I,  000. 

HCN 

30.  0 

200. 

HF 

30.  0 

100. 

so2 

20.  0 

500. 

H2s 

10.  0 

400. 

nh2 

500.  0 

2,  000. 

NOx 

25.  0 

200. 

coci2 

5.  0 

50. 

These  levels,  obtained  from  the  open  literature,  have 
been  selected  ag  estimates  only  for  the  purposes  of 
exercising  the  computer  simulation  using  these  levels 
and  are  not  intended  to  constitute  a conclusion  about 
the  actual  toxicological  effects  of  these  gases. 

To  further  compare  the  results  of  various  simulation  runs, 
a rough  estimate  of  the  combined  effect  of  the  toxic  gas 
concentrations  is  made  by  summing  the  ratios  of  each 
gas  concentrations  to  the  corresponding  irritation  or  life 
danger  levels  given  above.  If  the  sum  of  the  ratios  for  the 
irritation  level  exceeds  one,  the  indication  is  that  the 


40 


combined  effect  of  the  two  or  more  toxic  gases  may  cavsx 


short  exposure  irritation  even  though  no  single  concentration 
exceeds  its  irritation  level.  If  such  a condition  occurs,  a 
"1"  is  printed  to  the  right  of  the  statement,  "SE  IR  COMB 
EFF",  below  the  report  of  the  concentrations  for  each  gas 
distribution  condition.  An  analogous  computation  is  dot$ie 
for  the  life  danger  levels  and  is  reported  as  "SE  LD  COMB 
EFF".  Again,  it  must  be  emphasized  that  these  computations 
are  for  comparison  purposes  only  and  should  not  be  used 
for  any  other  purpose. 

3.  Flame  Spread  Data 

The  characteristics  of  each  distinct  fire  are  printed 
as  they  appeared  prior  to  the  flame  spread  calcula- 
tions for  this  time  interval.  The  remainder  of  the 
data  represents  conditions  as  they  existed  after  flame 
spread  calculations.  This  data  includes  the  number 
cf  flaming  and  smoldering  elements  on  each  surface, 
and  a summary  by  material  type  of  the  area  in  the 
flaming  and  smoldering  states.  In  addition,  a two 
dimensional  diagram  cf  each  surface  is  included 
which  enables  the  user  to  picture  the  state  of  each 
element  on  the  surfaces  involved. 

SAMPLE  OUTPUT:  The  following  five  (5)  pages 
contain  sample  output  from  a simulation  run  in 
which  IPEMS  = IP SPR  -10  seconds.  This  output 
resulted  from  the  sample  input  data  presented  in 
Section  4.  2 . 
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riHE=  40  SEC  AFTER  IGNITION 


FFF 


TIHE  = AO  SEC  AFTER  IGNITION 


SECTION  6 

PROGRAM  STATISTICS 


The  following  program  data  pertains  to  runs  made  on  the  Control 
Data  6600  Computer  System  installed  at  Area  B,  W right- Patterson  Air 
Force  Base,  Dayton,  Ohio  ( Building  676). 

Programming  Language: 

Operating  System: 

Computer  Storage  Required: 

Compile  Time  (CPU  seconds): 

Execution  Time  (for  sample  input 
run  to  7 30  sec.  simulated  time): 

Number  of  Cards,  Program  Source: 

Number  of  Cards,  Input: 


FORTRAN  IV 
NOS/BE 

120,000  words  (octal) 
18  seconds 
360  seconds 

approx.  3000 
approx.  150 
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SECTION  7 


* ’ ■“** 


PROGRAM  AVAILABILITY 


Because  of  the  length  of  the  code  of  the  DACFIR  program,  a listing 
of  the  code  has  not  been  included  in  this  volume.  Copies  of.  the  program 
code  and  sample  input  data  may  be  obtained  from  any  of  the  organizations 
listed  below.  * ' 


University  of  Dayton  Research  Institute 

Applied  Systems  Analysis  Section 

Attn:  Mr.  C.  D.  MacArthur 

300  College  Park 

Dayton,  Ohio  45469 

(513) 229-3921 

Department  of  Transportation 
Federal  Aviation  Administration 

Cabin  Fire  Safety  Research  & Development  Program 

Attn:  Mr.  C.  C.  Troha  ARD  520 

Trans  Point  Building 

2100  Second  Street,  S.  W.  , Room  1400 

Washington,  D.  C.  20591 

(202) 426-8416 

Department  of  Transportation 
Federal  Aviation  Administration 
Aeronautical  Center 
P.  O.  Box  25082 

Attn:  Mr.  James  Gillespie  AAC210 
Oklahoma  City,  Oklahoma  73125 
(405) 686-4374 


